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I show that a spontaneous electric polarization exists in the solution of the Falicov-Kimball model by mapping 
the strong coupling limit of this Hamiltonian into an xxz spin 1 /2 model with a magnetic field. In this way, I 
determine the phase diagram of the strongly interacting model and show the existence a transition to a mixed 
valence regime containing two phases: an orbitally ordered state and a Bose-Einstein condensation of excitons 
with a built-in electric polarization. 

PACS numbers: 71.27.+a, 71.28.+d, 77.80.-e 



The Falicov-Kimball Ql model (FKM) was introduced to 
explain semiconductor-metal transitions and has been exten- 
sively used to describe valence transitions in heavy fermion 
compounds. Its original version contains a dispersive band 
of itinerant d electrons interacting with localized / orbitals 
via an on-site Coulomb interaction. If hybridization between 
both bands is included, the / charge occupation is no longer 
a good quantum number, and it is possible to build coherence 
between the d electrons and the / holes. Based on a mean field 
solution of the FKM with a hybridization term, Portengen et 
al 0] proposed that this coherence gives rise to a spontaneous 
electric polarization associated with a Bose-Einstein conden- 
sate (BEC) of d — / excitons. 

Ferroelectrics are of considerable theoretical and techno- 
logical interest because of their highly unusual properties 
The ferroelectric (FE) transitions have traditionally been con- 
sidered as a subgroup of the structural phase transitions. As in 
the case of superconductivity, the existence of ferroelectrics 
based on a purely electronic mechanism would provide a new 
set of physical properties |2| and technological applications; 
for instance, it would open the possibility of controlling opti- 
cal properties with magnetic fields. 

The proposal of electronic ferroelectricity in the FKM was 
recently tested theoretically using different techniques. An 
analytical calculation in infinite dimension for the weak cou- 
pling limit lljj] did not confirm its existence. By using nu- 
merical methods to solve finite-size chains, Farkasovsky 0] 
arrived at the same conclusion for the intermediate and the 
strong coupling regimes. Recently, Zlatic et al 0] calculated 
the spontaneous polarization in the FKM from its exact so- 
lution in infinite dimensions. They found that the sponta- 
neous hybridization susceptibility diverges when the temper- 
ature goes to zero, indicating a possible non-zero polarization 
of the ground state. 

Hybridization between the bands however is not the only 
way to develop d — f coherence. An f — f hopping also 
induces it. Furthermore, I will show that in the mixed va- 
lence regime, the ground state of the FKM with / — / hop- 
ping is either an orbitally-ordered (chessboard) state or a BEC 
of electron-hole pairs (excitons). In particular, the condensate 
has a built-in spontaneous FE or antiferroelectric (AFE) po- 
larization induced by a pure electronic mechanism. The FE or 
AFE character of the solution depends on the relative sign of 
the d — d and the f — f hoppings. 

I obtain the complete phase diagram of the extended FKM 



by mapping the original Hamiltonian into an effective spin 
model. The mapping is exact in the strong coupling limit and 
the resulting spin Hamiltonian is a spin 1/2 xxz model with 
an applied magnetic field along the z axis. This spin model is 
exactly solvable in one dimension and its phase diagram has 
been determined very accurately for two dimensional {D ~ 2) 
systems (the exact solution has been numerically obtained for 
a 96 X 96 square lattice 01)- In this way, the results obtained 
in this paper prove that the phase diagram of the extended 
FKM contains FE and AFE phases induced by an electronic 
mechanism. 

Recently, high dielectric constants were observed in ox- 
ides of the type ACu3Ti40i2. In particular, the largest di- 
electric constant ever observed is exhibited by CaCu3Ti40i2 
(en ~ 80000 for single-crystal samples at room temperature) 
. In addition, high resolution x-ray and neutron pow- 
der diffraction measurements of CaCu3Ti40i2 rule out a con- 
ventional FE structural phase transition. The electronic ferro- 
electricity proposed by Portengen et al Ql was also ruled out 
due to the large value of the optical gap (A ~ l.beV) fllll . 
The strong coupling theory introduced in this paper shows that 
the excitons condense in the presence of a large gap. 

I will consider an extended FKM for spinless fermions on 
a D-dimensional hypercubic lattice: 

H = ed^nf + efY,n( + td^dldj 

i i (iJ) 

i (iJ> 

where nf — dj d;, and n( = f^f^ are the occupation numbers 
of each orbital. For historical reasons, I denote the orbitals by 
/ and d but in general they can represent any pair of atomic 
orbitals with different parity. For instance, in CaCu3Ti40i2 
there are d states which are hybridized with p bands [ 1 1 J . This 
fact is another important motivation to include a non-zero tf. 

The Hamiltonian H can be rewritten as an asymmetric Hub- 
bard model if the orbital flavor is represented by a pseudospin 
variable. A spin 1/2 is required to describe the two orbitals 
(flavors) on each site: 

4i = fl (2) 
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where the pseudospin generators are: 

The expression for H in the new language is: 



(3) 



+ C//''^nitnu + B,^rf, (4) 

i i 

where = ^(cd + £/) and = ed — ej. The new version of 
H is a Hubbard model with different hoppings for each spin 
flavor, t| = td and t| = tf, plus a Zeeman coupling with 
a magnetic field B^- Both terms break the SU{2) symmetry 
of the original Hubbard model (t^ = and = 0). The 
remaining symmetries are the U{1) groups associated with 
the conservation of the total charge and the total r^. In the 
original language, these U{1) symmetries correspond to the 
conservation of the total number of particles in each band 



0). 



I will consider from now on the half filled case, i.e. one par- 
ticle per site. For this concentration, it is well known that the 
Hubbard model in the strong coupling limit can be reduced to 
an effective Heisenberg model. In a similar fashion, I can re- 
duce H to an effective spin model when *C U-^'^. The low- 
est energy subspace for infinite U^^ is the one generated by 
states having one particle at each site, i.e., the charge degrees 
of freedom are frozen (the system is a Mott insulator) and an 
effective spin is localized at each site. In this limit there is 
a complete spin degeneracy because the energy does not de- 
pend on the orientation of each spin. To lift this degeneracy it 
is necessary to consider the lowest order processes in tc/U^'^. 
This can be done by a canonical transformation which elimi- 
nates the linear terms in the hopping and keeps the terms of 
quadratic order. Up to an irrelevant constant C ~ —NZJz/8, 
where Z is the coordination number and N is the number of 
sites, the resulting effective spin Hamiltonian is lll2il : 



(ij> 



j ' 



(5) 



. 2(t;+t?) 4tjt, . -1,0 

with J, = — ijTT^ and J± = -jjjt- J^eff is a spm 1/2 xxz 
model with an applied magnetic field along the z direction. 
The model is Ising-like ( > J±). However it is important 
to consider the whole phase diagram because the ratio Jz/J± 
can take any value if non-zero nearest-neighbor repulsions are 
added to H. 

Hgf f has been exactly solved in one dimension by means 
of the Bethe anzatz technique fisll . The one dimensional 
(ID) quantum phase diagram is similar to the 2D one that 
I describe below. The only important difference is that, as 
required by the Mermin- Wagner theorem IT^ . the excitonic 



condensate is critical at zero temperature (power law correla- 
tions) for the ID case. The phase diagram of i?e// has been 
determined recently for 2D systems by solving up to 96 x 96 
lattices with quantum Monte Carlo loop algorithm R|. The 
zero temperature phase diagram is shown in Fig. 1 and the 
corresponding name of each phase translated back to the orig- 
inal language of the FKM. Since H^f / is symmetric under a 
reflection in the xy plane, the phase diagram must be sym- 
metric under a change of sign of the magnetic field. The 
fully polarized solutions, obtained for large values of |i?z|, 
correspond to a full / band for positive B^ {cf ^ Cd) and 
a full d band for B^ negative (e/ ^ e^). The spectrum 
of both phases has a finite charge transfer (pseudospin gap) 
Act = \Bz — Z{\J±] + Jz)/2|. This gap vanishes at the 
quantum critical points \B'^\ ~ (|^±| + Jz), which are the 
boundaries for the mixed valence phase that emerges when the 
/ and d bands are sufficiently close: |e/ — e^l < ^{\Ji_\ + Jz)- 
If Jz > J_L, two phases are possible within the mixed valence 
regime. For small values of \Bz\, the Jz term dominates and 
induces a longitudinal antiferromagnetic (AFM) phase (chess- 
board state in the original language). When |i?2 1 is larger than 
a critical value, the magnetic field suppresses the Ising-like 
ordering and the J± term induces a magnetic ordered state in 
the xy plane(BEC of electron-hole pairs). The line separat- 
ing the orbitally ordered state and the BEC corresponds to a 
first order transition. This line ends at the Heisenberg point 
(J^ = J±) where both phases coexist. For Jz < J±, the only 
phase in the mixed valence regime is the BEC. In a real ma- 
terial, changing Bz = £d — can be achieved by applying 
pressure or alloying. 
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FIG. 1 : Two dimensional quantum phase diagram of H^f / obtained 
from Ref. 0|. The small circle indicates the position of the Heisen- 
berg point. The dashed line denotes the quantum phase transition 
between the mixed valence (non-shadowed) and the non-mixed va- 
lence regimes. 

The d — f exciton condensate has a built-in electric polar- 
ization [2]. This can be easily seen by realizing that the spin 
version of the order parameter for the BEC is the uniform xy 
magnetization for negative J± (condensation at k = 0) and 
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the staggered xy magnetization for positive J± (condensation 
at the AFM wave vector k = Q): 

= ^(rf X + y) for J_i < 

M^T = ^e*'^-'(Tf x + rf y) for > 0. (6) 

i 

Since is a two dimensional vector, it can also be repre- 
sented by a complex number |M-'-|e*'^ (tan0 = M^/M^) 
which is the usual expression for the order parameter of the 
BEC. The sign of J± is determined by the relative sign of td 
and tf. On the other hand, the uniform polarization operator 
is H: 

p = (^E(4/i+rfi//) = f (7) 

i 

where fi is the inter-band dipole matrix element and is the 
volume of the system. Therefore, the condensate has a built-in 
electric polarization which is proportional to the x (real) com- 
ponent of its order parameter For positive J±, the condensate 
becomes AFE because the staggered electric polarization: 

P5T = g E -'""-''(dlfi + djl) = ''-^M^sT (8) 

i 

is proportional to Mgj,. 

The three dimensional (3D) quantum phase diagram of 
Hef f fl4ll is similar to the two dimensional (2D) one shown 
in Fig.l. The same is not true for the finite temperature phase 
diagrams due to the Mermin- Wagner theorem llTsIl . The tran- 
sition temperature associated to the BEC is finite only for 
the 3D case. From the finite temperature phase diagram ob- 
tained in Ref. (tI, the BEC of electron-hole pairs undergoes a 
Kosterlitz-Thouless phase transition in a 2D system. 

I will now analyze the effect that a time dependent electric 
field Ee"^* induces in our FE or AFE condensate of excitons. 
The coupling term between the electric field and the uniform 
polarization: 

Hj=E-f,—Y,{dlf,+dJ^), (9) 

i 

corresponds, in the spin language, to the application 
of a uniform time dependent magnetic field Bi(t) = 
2E ■ fj, e*'^*x/51. From the point of view of the spin vari- 
ables this is like a magnetic resonance experiment since H^f / 
already includes a uniform static field B^z. Therefore, the 
equivalent magnetic system will have a resonant absorption 
at the frequency which tends to ujq ~ B^/h when Bi tends 
to zero. Back to the original language, this means that for 
small electric fields the optical absorption will be resonant at 
^0 = — £/ ■ This is an experimental fingerprint of the 
excitonic condensate. 

Since the above theory is only valid in the strong coupling 
limit, it is natural to ask whether the chessboard ordering 
and the excitonic condensate survive in the intermediate and 



weak coupling regimes. To answer this question it is more 
convenient to use the Hubbard-like representation of H (see 
Eq. 0). The dispersion relation for the non-interacting part 
of H {U^-^ = 0) is: e(k, (t) = + B,n + 2t, Y.u cos{K). 
At half filling and for = the Fermi surface of the non- 
interacting problem nests at k = Q. This indicates that an in- 
finitesimal value of U^'^ is sufficient to induce an AFM (chess- 
board ordering in the original language) instability. Again the 
presence of a non-zero magnetic field will induce a transition 
from the orbitally ordered state to the BEC of excitons. For the 
3D case, if C//'^ < J//'* - 2.85(|ia| + |ifc|) the magnetic field 
Bz induces an insulator-metal transition before the saturation 
of the magnetization (non-mixed valence regime) is reached 
fl^ . This means that for weak coupling a new metallic phase 
appears between the BEC of exctions and the non-mixed va- 
lence regime (see Fig. 1). For these reasons, I expect the phase 
diagram of H to contain an electric polarized BEC in the weak 
and intermediate coupling regimes as well. This subject will 
be studied more extensively in Ref. [17]. 

What happens if we consider electrons instead of spinless 
fermions? In this case each orbital can be occupied by two 
electrons and therefore it is natural to include local Coulomb 
repulsions U^^ and C/'*'^. By doing so the FKM is replaced 
by a two orbital Hubbard model and the large U"-'^ expan- 
sion gives rise to a Kugel-Khomskii like model fT^l contain- 
ing spin s (magnetic) and pseudospin t (orbital) degrees of 
freedom. If ta = -h = t and U^f = = U (the most 
general case will be analyzed in Ref. fl9il l the effective spin 
Hamiltonian is: 

Heff = ^Y.^i-S, + iJ,-Jo)Y,r.^r?+B^Y.^.: 

(iJ> (ij> i 

+ 2j2[Jor^r^ - Mrtr? + rfrj^)](Si ■ Sj + 1), 

(iJ> 

where Jo = At^ /U . The ground state of ^^e// ferromag- 
netic (FM) if Jo < Jq . For the FM solution, H^'^^ reduces to 
Hef f because all the electrons have the same spin orientation 
and hence can be considered as spinless fermions. Therefore 
the charge degrees of freedom (t) of the FM solution are ex- 
actly described by H^ff, and the phase diagram is the one of 
Fig. 1; i.e., ferromagnetism coexists with chessboard ordering 
or a FE BEC of excitons. If Jq > Jg, the system becomes 
AFM. In this case the effective transverse coupling for the 
pseudospin variables, J^'^-' = —2 Jj^((Si-Sj + i)), changes its 
sign because ((Si ■ Sj + |)) turns to be negative. For this rea- 
son the AFM solution coexists with an AFE condensate. This 
means that when a magnetic field induces a transition from 
an AFM phase to a FM one, it simultaneously changes the 
electric polarization from AFE to FE. This important property 
gives rise to new technological applications. 

From now on I will continue with the spinless case just to 
isolate the basic mechanism for FE and AFE which is associ- 
ated with the charge degrees of freedom. An important aspect 
of this analysis is the inclusion of a non-zero hybridization 
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term: 

Hy = VY,{dlf,^,^+fl,J,-f^d,^-_^-dl,J,), (10) 
1,1^ 

where i/ runs over the different spatial directions (x, y and z 
in three dimensions). The different signs in the hybridization 
terms are due to the different parities of the two orbitals. (The 
crystal has inversion symmetry.) By adding Hy to H, we get 
the following additional terms for Heff (large U expansion): 

where J' W^/U^'^ and 4^ = W{ta + h)/Uf'^. Let 
me now consider the perturbative effects of the hybridization 
(J', ^ Jz, J±) on the phase diagram of Fig. 1 for the 
FE case (J± < 0). The first term is a Heisenberg interac- 
tion which just produces a renormalization of Jz and J±. The 
mean value of the second term is zero in the FE phase (BEC 
of excitons) and therefore they do not make any contribution. 
The last term introduces an easy axis anisotropy along the x 
direction and lifts the U(l) degeneracy. The BEC then is re- 
placed by an Ising-like FE state characterized by the breaking 
of the remaining Z2 symmetry. Therefore, the spontaneous 
ferroelecricity remains when the hybridization is included per- 
turbatively; however, the resonant response to a time depen- 
dent electric field (see Eq. ||9l) disappears due to the absence 
of Goldstone modes. In other words, the hybridization makes 
the electronically induced FE phase similar to the ones in- 
duced by structural phase transitions. 

In summary, I derived the phase diagram of the FKM with 
atf hopping term in the strong coupling limit. The insulating 



phase obtained at half filling has a transition from a non-mixed 
valence to a mixed valence regime as a function of the energy 
difference between the centers of both bands. Two different 
phases are present in the mixed valence regime: a BEC of ex- 
citons with a built-in electrical polarization which starts just 
at the valence transition and an orbitally ordered (chessboard) 
state which appears when the centers of the bands are suffi- 
ciently close. These results were extended to the intermediate 
and weak coupling regimes due to the nesting property of the 
Fermi surface of hypercubic lattices. I also mentioned the ef- 
fect of including spin degrees of freedom in H: the interplay 
between magnetic and charge degrees of freedom gives rise 
to the coexistence of FE and FM phases which are coupled 
to each other This opens the possibility of controlling optical 
properties by applying magnetic fields or controlling magnetic 
properties by applying electric fields. 

The effect of a non-zero hybridization was also considered. 
The main conclusion is that the U{1) degeneracy associated to 
the BEC of excitons is lifted by the hybridization and replaced 
by an Ising-like FE state (broken Z2 symmetry). Then, the 
resonant response to a time dependent electric field disappears 
because the Goldstone modes acquire a finite mass (gap). 

These results indicate that the following characteristics are 
favorable to the formation of an electronically driven FE state: 
a) The system must be in a mixed-valence regime and the two 
bands involved must have different parity, b) It is best, though 
not necessary, if both bands have similar bandwidths. c) A 
local Coulomb repulsion {U^'^) between the different orbitals 
is required, d) The hybridization between the bands must be 
small compared to their bandwidths. 
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